arXiv:1508.03224v3 [math.PR] 30 Jan 2016 


SOME PROPERTIES OF PRABHAKAR-TYPE 
FRACTIONAL CALCULUS OPERATORS 


Federico Polito, Zivorad Tomovski 


Abstract. 


In this paper we study some properties of the Prabhakar integrals and derivatives and of 
some of their extensions such as the regularized Prabhakar derivative or the Hilfer-Prabhakar 
derivative. Some Opial- and Hardy-type inequalities are derived. In the last section we point out 
on some relationships with probability theory. 


1. Introduction and background 


The aim of this note is to study the properties of some integral and differential 
operators that can be related to a specific convolution-type integral operator, called 
Prabhakar integral, introduced by Prabhakar ll24ll . Before describing the results obtained 
we start by recalling the basic definitions and mathematical tools that will be useful in 
the following. 

First, let us give to the reader some insights on the classical operators related to 
fractional calculus. For more in-depth information, applications and related topics it is 
possible to consult some of the classical references, e.g. Samko et al. Il25l . Podlubny 
12^ . Kilbas et al. ifTSl . Diethelm 18], Gorenflo et al. ifTSl . Mainardi We also 
suggest the reader to consult the references listed in these books as the relevant literature 
is becoming richer and richer. 

In the following we give the definitions of the Riemann-Liouville integral which 
in some sense generalizes the classical multiple integral, and of its naturally associated 
differential operator, i.e. the Riemann-Liouville derivative. We then proceed by describ¬ 
ing the regularized version of the Riemann-Liouville derivative, the so-called Caputo or 
Caputo-Dzrbasjan derivative, introduced independently in the sixties by Caputo 0121 
and Dzrbasjan and Nersesjan Qo). 

Definition 1 . (Riemann-Liouville integral) Let / — oo < a < t < 

b <oo^ he a locally integrable real-valued function. Let us further define the power-law 
kernel -JCait) = f“^Vr(o:), a > 0. The operator 



a > 0 


( 1 ) 


is called Riemann-Liouville integral of order a. 
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Definition 2. (Riemann-Liouville derivative) Let / £ L^{a,b), <a<t< 

b <o°, and S m = [a], a > 0, where W'"’^{a,b) is the Sobolev 

space defined as 


W'”’\a,b) = \feL\a,b):—fG L\a,b) \ . 


( 2 ) 


The Riemann-Liouville derivative of order a is defined as 

jm 1 d"* 

{D^^f)it) = -/^“/(O = ^^(3) 

V a+J ly ) df™ “ r{m-a)dt"'Ja 

We denote by AC" {a,b), n £ N, the space of real-valued functions /(f) with 
continuous derivatives up to order n — 1 on {a,b) such that /("^^) (f) belongs to the 
space of absolutely continuous functions AC {a,b), that is, 

AC"(a,h) = |/:(a,fi)^K: ^/(x) £ AC(a,h)|. (4) 


Definition 3. (Caputo derivative or regularized Riemann-Liouville derivative) 
Let the parameter a > 0, m = [a], and / £ AC'"{a,b). The Caputo derivative (also 
known as regularized Riemann-Liouville derivative) of order a > 0 is defined as 






^ ft 

-r / (f-s)'"-l-“—/(S)di. 

-a)Ja 


T(m — a) 


(5) 


The above Definition should be compared with the non-regularized case of 
Definition |2l The reader should also be aware of the fact that the above derivatives 
can be defined in different ways. For example a more interesting and intuitive way of 
defining the Caputo derivative is given by the following theorem. 

Theorem 1.1. For f C AC" {a, b), m= [a], a£K+\N, the Riemann-Liouville 
derivative of order a of f exists almost everywhere and it can be written in terms of 
Caputo derivative as 


Theorem 1 1.1 1 is interesting in that it in practice describes the set of functions with 
which the Riemann-Liouville derivative can be regularized. It is well-known that if 
/(f)£AC-(fl,h), 


(0=0, Q<k<m-1. (7) 

By taking the Laplace transform of both sides of (|^ we immediately realize that (|^ still 
holds if (0 is true. 
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In the recent years new alternative definitions of fractional operators have been 
introduced in the literature. An interesting example is the so-called Hilfer derivative 
OS Ell El. The idea behind the introduction of this derivative is to interpolate between 
the Riemann-Liouville and the Caputo derivatives. As it is clear from the definition 
below, the Hilfer derivative depends on the parameter v G [0,1] that balances the 
individual contributions of the two fractional derivatives. 

Definition 4. (Hilfer derivative) Let /r G (0, 1), v G [0, 1], / e l' [a,b] , —OO <C 
a < t < b < oo^ /* € AC' [a,b] , where J^(f) = ■ The Hilfer 

derivative is defined as 

Notice that Hilfer derivatives coincide with Riemann-Liouville derivatives for 
V = 0 and with Caputo derivatives for v = 1. 

In order to proceed with the description of the relevant operators involved, we need 
to consider now the function 


ep.M.o(0 f GM,p,M,®,7eC,91(p)>0, 


(9) 


where 


£■7 / N ^ y E{Y+k) x'^ 

t,r{r)ripk+^)k\- 


( 10 ) 


is the generalized Mittag-Leffier function first investigated by Prabhakar ll24ll . The so- 
called Prabhakar integral is constructed in a similar way of Riemann-Liouville integrals. 
The main difference is that, the power-law kernel in the integral representation of 

the operator is replaced by the function (|^. The kernel (|^ actually generalizes 
in the sense that gp p ^(f) = ■ The Prabhakar integral is hence defined as follows 

Giiini. 

Definition 5. (Prabhakar integral) Let f G {a, b), 0<a<t<b<o°. The 
Prabhakar integral is defined as 

{t -y)P]fiy)dy, (11) 


where p,p,a),7G C, with 91(p),91(p) > 0. 

As just remarked, for 7 = 0, the integral GD coincides with Q- 
An interesting property of the Prabhakar integral is the following ifTTl formula 
( 2 . 21 )]. 

Proposition 1.1. Let 7,p,p,ij,(7,a) gC, 91(p),91(p),91(Tj) >0, f gK. Then 

(Ep,p,„,or<r,.»)(0 = (12) 


Analogously to the classical fractional operators, a related differential operator can 
be defined as follows. 
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Definition 6. (Prabhakar derivative) Let / G L*(a,Z7), 0 < a < f < < 0°, and 
(b(’) £ m= [jLi]. The Prabhakar derivative is defined as 




/)(0 



(Ep 


-r 

,m—11,(0, 


/)) ( 0 , 


(13) 


where ii,co,Y,p gC, 9^(^),9^(p) >0. 

The inverse operator ( [T3| ) of the Prabhakar integral, for 7=0, generalizes the 
Riemann-Liouville derivative as Ep [(o{t —7)^] = 1 /r{m — p) in the kernel of ( [T3| l. 
See also Kilbas et al. ifTTl for a different but equivalent definition. 

The analogous operator to the Caputo derivative, that is the regularized Prabhakar 
derivative plays an important role in the construction of meaningful initial-value prob¬ 
lems as it was noted in Garra et al. 112. The regularized Prabhakar derivative was 
introduced in D’Ovidio and Polito 

Definition 7. (Regularized Prabhakar derivative) Consider / G AC"'{a,b), 0 < 
a < t < b < °o. The regularized Prabhakar derivative reads 

m—1 

k=0 

For 7 = 0 the operator ( [l4| ) coincides with the Caputo derivative (|^. 

Remark 1. Let p >0 and / G AC"'{a,b ), 0<a<t<b<oo, Then 

where h{t) = f{t) -1,'"-^} (a+). 

The last operator we need to derive the results described in Section]^ is similar 
to the Hilfer derivative but based on Prabhakar operators. Therefore, let us give the 
following 

Definition 8. (Hilfer-Prabhakar derivative) Let p G (0,1), v G [0,1], and 
let /G L'(a,fi), 0<a <t <b<oo, /* G . The Hilfer- 

Prabhakar derivative is defined by 




/)(0 


(A 

P.v(l-p).ffl.fl+df 


(E 


-r{i-v) 

p,{l-v){].-fi),(0,a 


+ 


/)) ( 0 , 


(16) 


where 7, ® G K, p > 0, and where (E° g ^ ^+/)(f) = /(f). 

This Hilfer-Prabhakar derivative interpolates the two Prabhakar-type operators ( fTB] ) 
and (|T^ and it specializes to the Hilfer derivative for 7 = 0. 
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We can also think of a regularized version of ([T^, that is, for / G AC* {a,b), we 

have 


= (E;kp«4-^) (0^ 

Notice that in the regularized Hilfer-Prabhakar derivative ( [T7| ) there is no dependence 
on the interpolating parameter v. 

Before proceeding to the results section we present here an estimate for the function 
/3 ft) proves to be necessary in the proof of some of the theorems below. We 

refer in particular to Theorem 3 of Tomovski et al. EH- 

Theorem 1.2. (Tomovski-Pogany-Srivastava) For all a G (0,1), 7,cu > 0, 
ay > p — I > 0, the following uniform bound holds true: 


'a.p,m 


(0 


< 


r(7-g^)r(g^) 


7ra®(/*^*)/“r(7) (cos(7ra/2))*' **/“’ 


t > 0. 


(18) 


The proof of Theorem 1.2 makes use of an estimate of the Wright function given 
by Stankovic 


2. Some bounds and operational calculus with Prabhakar-type operators 


The hrst two results we present concern the boundedness of the Prabhakar integral 
0 in the space , p G (0,1] and in L* for functions, p G (1,°°). 

Theorem 2.1. Let a g (0,1), 7 ,® > 0, and aY> p - 1 > 0. If (p g LP {a,b), 
0 < p < I, then the integral operator p m a+ bounded in {a^b) and 




<M\\(p\, 


(19) 


where the constant M, 0 < M < 00 , /s given by 


M = 




naco^a' [cos [^)Y “ 




in which Be(/i, v) is the Beta function. 

Proof In order to prove the result it is sufficient to show that 


( 20 ) 


(Ea./3.®.a+‘P)|r= / / 'Elp[o}{x-tY](p(t)At 

^ P J a J a ^ 


ck<oo. (21) 
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This can be done by recalling the well-known integral inequality 


I f{t)dt <[ |/(f)l'’df, 0</?<l, 

J a J a 


( 22 ) 


and the uniform bound of the function eL o ,At) (see Theorem 3 of Tomovski et al. 
ED). We obtain 


(23) 


(E’' « +(p) 

a.p,a),a+^ > 

"iiAr 

P Ja \Ja 



V 

(r-«7)r( 


^;raa)^a'r(7) [cos(^)]’' “ 
< [b — a] 


1^ (f)l^df ) dr 




nao) a* [cos (^)]^ “ 


l<Pll^ 


This completes the proof. 


□ 


Theorem 2.2. Let a e (0,1), 7,® > 0, and a'Y> j5 - 1 > 0. If (p e Lf {a,b), 
p > 1, then the integral operator ^ ^ is bounded in L} {a,b) cind 


(E’'fi +<P) 

^ a.p,)o,a+^-' 


<M 


(b — a) 


1 


q-\- 1 


1 /^ 


1 <PIL: 


(24) 


where \/p-\-l/q= 1. 


Proof. By Fubini’s theorem and Holder inequality, we obtain 


(E 


a,(3,CD,a 


(P) 


f I 


A-tY 


a,p 


[a)(x-f)“] ^(f)df 


dr 



(25) 


where the constant M is given by @. □ 

We now give a result on the boundedness of the Hilfer-Prabhakar derivative ( [T6] l in 
the space . 
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Theorem 2.3. For /i G (0,1), V G [0,1] and f G the operator 


is bounded in the space L} {a,b) and 


where 


Ml = [b-aY^^^^'^ ^ 


l(y(v-i))J 


i^^Qfr{pk + v{l-p))\lpk + v{l-p)] 

l(7(v-l))J 


> 0 , 


M2 = (b — aY^ ^ Y' -—-- 

^^^\T{pk + pv-p-v)\[pk + pv-p-v] 


> 0 , 


(26) 

(27) 

(28) 


and 7,0 S K, p > 0. 

Proof. Using the estimate given in Theorem 4 of Kilbas et al. ini, we obtain 

d 




p-rv 


.o.«+ df 


(29) 


<Mi 


= Mi 


d / 

Yt V®^p,(l-v)(l-Ai) 


-y(l-v) 


(0 

<MiM 2 ll/ll 1, 

where Mi and M 2 are the constants defined by (|27| and ( |2^ . □ 

We now proceed to the study of the boundedness property of the regularized 
Prabhakar derivative ([T^ and of the regularized Hilfer-Prabhakar derivative ([T7|i. In 
particular we prove L} boundedness. 

Theorem 2.4. If feW'"’^ (a,h), 7,® e M, m = [p], p > 0, then regularized 
Prabhakar derivative is bounded in l} {a,b) and the following inequality holds true: 


where 


_ /7 \in —p % ’ I( y) t; I 

^^^'\r{pk + m-W[pk + m-p] 


kl 


(30) 


(31) 


Proof Using the estimate for the Prabhakar integral operator (see Kilbas et al. 
ini), we obtain 


(Ct)! f) 


r—/ 


<K 


f 


(m) 


(32) 

□ 
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Theorem 2.5. if f {a,b), 7,® e M, /z e (0,1), p ^ 0. thcfi the regular- 
ized version of the Hilfer—Prabhakar derivative ^ is bounded in {a,b) and 

the following inequality holds: 



(33) 


where 


\i-r)k 


^ £|r(p^+i-M)|[p^+i-M] 


® (P —a)' 


k\ 


(34) 


Proof As in the proof of the preceding theorem we nse again the L} estimate for 
the Prabhakar integral operator ini, obtaining 




E 


-r _ 


t:/ 




(35) 

□ 


The following Proposition |2. l| and Theorem |2.2| present basic composition relations 
involving Prabhakar integrals and Hilfer-Prabhakar derivatives. Some specihe cases are 
highlighted. 

Proposition 2.1. The following relationship holds true for any Lebesgue inte¬ 
grate function (jd € L* {a,b): 




(36) 


where 7 ,5, ® G K, p,A> 0 , pG(0,l), ve[0,l], A ^ p -j- V — pv. Ill pQ.rticulQ.r, 

(^;;^’>(e^,a..«./)) = (4V"/)- (37) 

Proof Using the semi-group property of the Prabhakar integral operator ifTTl . we 
obtain 




( — 
P,v(l-Ai).ffl.a+df 




-r(i-v) 

p,(l-V)(l-p),(B,fl+ 



/p-yv ^fp-r(l-v)+5 

1 ^p,v(l-p),CD,a+ \^^p,(l-v)(l-p)+A,ffl,a 




-yv 

p,V(l-/l),C0,t3+ 


y p,(l-v)(l-p)+A-l,(H,a+ 


-/))m 

/)) W 


(E 


S-y 

'p,A-p,m,a+ 




□ 
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Proposition 2.2. The following composition relationship holds true for any 
Lebesgue integrable function ^ {ti,b): 


where 7 ,® S K, p,X >0, p G (0,1), V e [0,1], A > p + v-pv. 


(39) 


Proof It is sufficient to prove the first relation. The proof of the second follows the 
same lines. We have 




= E 


-yv ^ 

= ( — 
P.v{l-p),(B.a+ dr 


7 (^ v) jX „ 

V^P,(l-v)(l-p),( 0 ,a+^a+ 




r(i-v) 

'p,(l-v)(l-p)+A,a),a+ 


it) 

it) 


(40) 


(®'p,v(l-p),fl),a+ (®^p,(l-v)(l-p)+A-l,( 0 ,a+*(^) ) i^) 

= i^pl-,,a,,a-tp){t). 


On the other hand. 


(') 


= /; 


E-rv 

p.v(l-p),(B,a+ dr 


= ( — 
p,v(l-p)+A,CD.a+ dr 


V®p.(l-v)(l-p),ffl,a+^ 
V%.(l-v)(l-p)+A,ffl,a+‘P 


it) (41) 


(0 


(®p,v(l-p),(B,a+ (®^p,(l-v)(l-p)+X-l,(B,a+^ 

-y 

p,?i—fi,co,a 


= i^A-unr.^tp){t). 


□ 


Example 1. Asa didactic example we calculate the (non regularized) Hilfer- 
Prabhakar derivative of the power function r'’^*, p > 1, with a = 0. As in Definition|^ 
we consider p G (0,1), v € [0,1], 7 ,® G K, p > 0. We obtain 




(42) 


= E 


-yv 


E 


-y(l-v) 


'p,v(l-p).CD. 0 +dr V P,(l-V)(l-p),(B,0+ 


tp- 


(x) 


=^ip){^p^p 


A (Ai-v)ii-p)+p-i^-r{i-v) (cotP)]\ (x) 
-p).ffl,o+drV %,(i-v)(i-p)+p(®f 


= r 
= r 


E 


j(i-v)(i-p)+p-2E-r(i-v) 


-yv 

'p,V(l-p),CD,0+‘ ^p,(l-v)(l-p)+p-l 

v(i-p)-iE-yv 


(x-r)''^^ P> 'EpJ^^^_^^icoix-t)P) 

= xP-P-^EpJ,_^{coxP). 
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In the last step we applied Proposition [TTT] 

Example 2. Similarly to Examplewe calculate the (non regularized) Hilfer- 
Prabhakar derivative of the function ^ ^ (f). We set € (0,1), V G [0,1], 7,0) G K, 
P >0, j3 > 1. In this case we have 




(43) 


— ( 2!_ 


- ( 

- ' *^p,v(i-p) 


d 

df 

d 


,®,o+^p,/t,ffl (o)) w 


A /pi-v)(i-p)+i3-i£yv (mfP)^ 1 (x) 

,»,o+df V %.(i-v)(i-p)+i3 I) J W 

Jl-v)(l-p)+i3-2prv (mtP)\(x) 

— r/^-p-lpO ('( 3 „P'i _ J*_ 

As in the above example in the second-to-last step we applied Proposition |1.1| 


3. Opial- and Hardy-type inequalities 


Opial-type inequalities have a great interest in mathematics in general and in 
specific fields such as the theory of differential equations, the theory of probability, 
approximations, and many others. The classical Opial inequality was introduced by 
Opial fT2\ and reads as follows: 

Proposition 3.1. Let f{y) g C'(0,/r), /(O) = f{h) = 0, and f{y) > 0, y G 
(0, h). It holds 

1 /( 7 )/'(7)1 dy < J I’' {f'iyfdy, (44) 


where hjA is the best possible constant. 

During the past years the classical Opial inequality has been generalized by many 
authors in many different directions. Eor a thorough account the reader is suggested 
to consult the monograph on Opial-type inequality by Agarwal and Pang m. See also 
the recent paper by Earid et al. ca. Eor further Opial-type inequalities involving the 
classical fractional operators see Anastassiou la¬ 
in this section we describe and prove several Opial-type inequalities involving as 
differential and integral operators those analyzed in the previous section. 

Theorem 3.1. let f e l\0,x) . If a e (0,1), 7 , m > 0, a 7 > j3 -1 > 0, p,q> 
1, l/p+l/q = I, then the following inequality holds true: 



y 

a,p,afi+ 


f) 


x^Ip 

(s) \f{s)\ds<K— 



(45) 
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where 


Be (7-^,^ 


Ttaco'^[cos “ 

Proof. Applying Holder inequality to the Prabhakar integral operator, we get 




, \ r 

I (m) du 


\f{u)\'^du 


r(7-^)r(2^) / « X 1/9 

< - ^ ^ J \ ^ J f/p( / |/(M)|^d«) . 

7taco^r{Y)[cos(^)Y^^ ^ ^ 

Let z(s) = /o |/(M)|‘^dM. Then z'(s) = |/(s)|^, i.e. |/(s)| = Hence, 

Applying again Holder inequality, we obtain 

I |K,/3,«,or/)(")|l/(")|d*^^Q n {z{s)f{s))ds\ 


_ x^/p (zix)Y^‘> _ X^P f F 

21 /p 21/9 2 Wo 


|/(s)|'idi 


The following theorem concerns an Opial-type inequality involving at the same 
time both Hilfer and Riemann-Liouville derivatives. 

Theorem 3.2. Let JJ. € (0,1), V € (0,1], f & L} (0,x), x > 0. Let furthermore 
f) € L°° (0,x). Then, for 0 < p < 1, the following inequality holds: 


Wf) If) (I) I * > ©(-v) (f\[ 




where 1/ p+l/q = 1 and 


0 (x) = 


{v{l-n)p-p+2)/p 


r(v(l-Ai))((v(l-Ai)-l)p + l)'^'’ (v(l-p)p-p + 2y/P 


Proof Recall that the Hilfer derivative can be written as 


K;7)(0 = (4/'‘’"^''7)) f) = (^^“'(0^2"-'"/)) (,) (48) 
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Applying now the reverse Holder inequality (remember that here p € (0,1)), we obtain 


IK^/) (01 


> 


1 


r(v(i-At)) \Jo 
1 


t \ i/p 


f({v(l-Ai)-l)p+l)//) 


r(v(l-Ai)) ((v(l-^)-l)p + l)i/'’ \Jo 


D: 


I Do. 

iM+v-AiV 

0 + 


Ll + V-jlV 


f) (t) ^dTj 

l/q 


(49) 
\ i/<? 


dr 


Let us 


now define z(f) =/q (<I>(s))'^di, where <I>(f) = (0 


4>(f) = (z'(f))'/^i.e. (0| =(0(0)'/'- Hence, 




0 + 


£»: 


0 + 

0 + 


dt 


> 


1 


(z(f)0(f))'^'^df 


> 


rx f{{v{l-^l)-l)p+l)/p 

r(v(l-Ai))io ((v(l-^)-l)p + l)*/P 

1 

r(v(l-Ai))((v(l-Ai)-l);^ + l)'/'’ 

)iMi-p)-i)p+i)/py^f'^ ^''^''(z(f)z'(0)df^ 

1 ^((v(l-|j)-l)p+2)/p 

r(v(l-Ai))((v(l-Ai)-l)p+l)‘/'’ {(v(i-m)-i)p+2)‘/'’ 21/? 

2-1/9 


r(v(i-Ai))((v(i-Ai)-i);5 + i)‘^'’ 

jf((v(l-tJ)-l)p+2)/p 


((v(l-p)-l)p+2)>/'’ IJo 


d: 


,p+v-pv 

0 + 


\ 

df , 


Then 


(50) 


and this proves the claimed formula. □ 

Corollary 1. Let/re (0,1), /eL^(0,fi), feL{0,b), 0<x<b. For 
0 < p < 1 the following inequality holds true: 


£ I {^D^f) it) • /' (f) I df > 0 (x) I/' (f) I df) , 


(51) 


where 


0(x) = 


2 -“? 


r(i-/i)((-/ip + i)(-/ip+2))'/'’ 


(-pp+2)/p 


ant/ ^ ^ = 1 • 
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Theorem 3.3. Let f G {0,x ), x>0 , f G AC™ (0,.*) and (D^ ^ ^ q+/) e 

L°° {0,x), 0 < p < 1, 0) > 0, 7< 0 , (0+) =0, ^ = 0,1,2,... ,m — 1, anrf —P7> 

m — — p,^>l, 1/p+ 1/^=1 f/ien the following inequality holds true: 


[ (') (d;,„,„,o./) (>) I d< < Si («) (//1/<"> (.) 


where 


LI (x) = 


r(- 7 - 


m—/i —1 1 p / m—/i —1 




r2/p 


Ijq 


(52) 


m—U—1 ,/ m—U-1 r\ 

TipCO P r ( — 7) [cos (^)] ^ 


Proof. Using Holder inequality and the uniform estimate of the function e, (see 
Theorem 3 of Tomovski et al. ED), 



where 


M = 


m-u-1 m-u-l 

Tzpco p r(-7) [cos (^)] ^ P 


Let z{t) = Jq (i) di. Then z'(f) = f^"'^{t) , i.e. (f) = (z'Hence, 


/ 


/W (f)(D^ 




/) (f)|df (z(t)z'(t))‘/^df 




i/g 


= M 


x^Ip {z(x)fl‘i x^/P 


2i/p 2 '/? 




\ 2/? 


□ 


The following two results are given without proof. However, the claimed formulae 
can be derived with similar methods to those implemented in the proofs of the previous 
theorems of this section. 
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Corollary 2. Let / C L' (0, y) , y > 0. Let furthermore f G AC”^ (0, a) , m C N, 
^ 0 +/) ^ L’^iO^x), 0 < p <1, CO >0, 7<0 and -pY> m-p-l > 0. If 
p,q > I, l/p+l/q — I, then the following inequality holds true: 





where 


Q. (a) 



m—n — l 
P 


m—^ —1 
P 


m—U—l f., 

Tzpco P r(-7) [cos (^)] ^ P 


^2/p 

~Y~ 


(53) 


Theorem 3.4. Let / e L'(0 ,y), x > 0, ^ e AC'(0,fe). Lef 

furthermore (0^^^+/) G L°° (0,y), /i S (0,1), V £ [0,1] ,0 < p < 1, (0 > 0, 7 < 0, 
and —py > I — p > 0. If p,q > I, l/p+l/q = I, then the following inequality holds 
true: 


‘ p,(O.0+ 


df V 


< Q. (x) 


d /„-y(l-v) 


- fv-' 

df V P’(i-''){i-p),ffl,o+ 


df 


■7 \ 2/9 

dfj , 


(54) 


where 


Q, (x) = 


r(-7v- 


v(l-p) ) p { v(l-p) 


r2/P 


;rp®'^r(-7v)[cos(f)]-’'''-^ ^ 


> 0 . 


Theorem 3.5. (Hardy-type inequality) Let p,q>l, \/p+l/q=l, (X,j5,y,co > 
0. If f G L‘1 {a,b), a < b, then the following inequality holds true: 


where C = 


[\{K,p,<o,a^f) (ofdf<c/j/(f)rdf, 

G (0,1), a7> j3- 1, we have 


(55) 


/ I {K,p,<o,a^f) ^ (56) 


where K = M {b — . The constant M is given by the right hand side of ( |18| ). 

Proof We prove only inequality (|55|l as the proof of (|56ll follows the same lines as 
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the proof of Theorem |2.1| By applying Holder inequality we have 


l/(T)rdT 


/ rt P / r‘ \ 

<(/Kf,10^) ( 

^4/3+1,l/(0l^df 


< 


'(/ i/(ordf), 


Thus we have 

(eL/3,».„+/) (0 

for every f € [a,/?]. Consequently we obtain 


< 


df^ l/(0l'^df 

I 1/(0Tdf 


(57) 


(58) 


(59) 


'«,/3+2,fi,(^-«)]^/ i/(ordf. 

In the last step we made use of formula (5.5.19), page 100, in Gorenflo et al. ca. □ 


4. Applications to probability theory 


In this section we make some remarks about connections of the described operators 
with the theory of probability and of stochastic processes. 

First we should recall that some applications to probability have been already 
discussed in several articles. As an example we mention the paper by D’Ovidio and 
Polito 13 in which the regularized Prabhakar derivative has been defined with the purpose 
of introducing a broad class of stochastic processes related to some partial differential 
equations of parabolic and of hyperbolic type. The paper by Garra et al. ca instead, 
presents an analysis of a generalized Poisson process in which the governing difference- 
differential equations contain a regularized Prabhakar derivative. The interested reader 
is encouraged to consult also the references therein. 

Consider the Wright function (j){—a,p;z), a G (0,1), p,z G K, defined as the 
convergent series 


^{-a,p-,z) 


y _^_ 


This can be used to construct a probability density function which proves to be interesting 
in many aspects. We will refer basically to the papers by Stankovic 1291 . by Gorenflo 
et al. lT4l . and by Tomovski et al. 1311 . for the properties of Wright functions. 
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We know [see also [18] O that the Mellin transform of 0 for the following specific 
choice of the parameters reads 

[3-a7;-z)) (7) = a,j3 G (0, 1), 7> 0, j3 > a7. (60) 

From this, it is clear that we can consider a random variable, say X , supported on IR+, 
such that its probability density function is 

g{x) = -7a;-x)l(o,„,)(;c). (61) 

and such that EX = 7r(j3 )/r(a + )3). The above definition is justified by the positivity 
of ( |6 T] i as remarked in the proof of Theorem 2 in Tomovski et al. OTl . 

The above random variable is particularly interesting in that it generalizes the 
marginal law of an inverse stable subordinator. Let us thus consider an a-stable 
subordinator y“(f), f > 0, that is an increasing spectrally positive Levy process such that 
Eexp(— Ay“(f)) = exp(— fA“) Let us call £“(?)= inf{i >0: y“(i) ^ (0,f)}. 

The marginal probability density function f{x,t) = P(£’“(f) G dx)/dx satisfies the 
fractional pde 


(^D“+J){x,t) = -^f{x,t), x>0,t>0, (62) 

and can be explicitly written as 

/(;c,f) (-a,l-a;-^) , ;c> 0, f > 0. (63) 

The function ([6T]i can be rewritten by considering a fixed time t as 

(-a,i3-a7;-^) , x>Q,t>Q, (64) 
which clearly generalizes © and ( j^ . 

Theorem 4.1 . The space-Laplace transform of the probability density function 
© writes 

Moreover its space-time-Laplace transform reads 


ajccr-li 

|(s,t(7) 5>0, CI>0. (66) 
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Proof. By using the contour integral representation on the Hankel patlj^of the 
Wright function (j) we obtain 


g{s,t)= [ e?-F?“C“^-/^dCck (67) 

Jo r( 7 ) 2m JHa 


nii)r 


ya 


r(7) 2%i Jna 

mr 




ya 


eii^ay-p 


2ni 

m) I 


e^^ay-p 


.dC 


dC 


2m ina (f“i + C“)’' 

= r(^)£^^(-5f“). 

The last step is justihed by using the contour integral representation of the reciprocal of 
the Gamma function. 




r(rj) 


and by the following calculation: 


</3W = E 


zT( 7 +r) 


^^T{Y)r\T{ar + fi) 
1 ^zT(7+r) 


( 68 ) 


(69) 


J_ y z’T{y+r) r ^t;r-ar-p^r 
2m r( 7 )r! jHa 


2m jHa 


c-'^oc E 


/ 7 +r -1 


r =0 




(zC-y 


JRa. 


r =0 




1 

2m 


,f e?C^^l-2r“)^^dC 

I JRa 


1 r 

~ 2m Jh, (C“-z)’' 

The space-time-Laplace transform can now be easily obtained as 

djay-P 


|( 5 ,CT)=r(/ 3 ) 


(n7“ + i)r’ 


s > 0, d > 0. 


(70) 


□ 


^The Hankel path starts at (— 0 °, —£), e > 0, proceeds to the origin on the lower half-plane, circles the 
origin counterclockwise and then returns to (—oo,e) along the upper half-plane. 
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Theorem 4.2. The function g{s,t) satisfies the fractional equation 

, , fP-ay-l 

[{D^Tr+^Ys) (^>0 = ^^P'^ rip-arY " > ^ 

Proof The proof follows from the properties of the time-Laplace transform of 

I- □ 

Remark 2. Notice that, for j3 = 1, the function § and the governing equation 
•ED are connected to the Havriliak-Negami relaxation (see Stanislavsky et al. Il28l for 
more detailed information). 

Let us now recall Theorem 1 of Tomovski et al. ED 

Theorem. (Theorem 1, Tomovski et al. ED) For all a € (0,1], j3 > 0, 7 > 0, 
t > 0, we have 


1 p 

where T£t stands for the Laplace transform with Br^j^ the Bromwich path (i.e. {i = 
a + iT: a > (To, T G M}) and 


KlM = 


r^r-P (rarctan " “^0 


71 (r2“+2r“cos(7ra) + 1)^^^ 

Furthermore, for all a G (1,2], j3 >0, and 7 = n G N, 

/ ^an-fi \ 2(-l)"-' 


^a./3.l(0=^r 


1 


(s“ + l)”/ a"(n-l)! 


fcos(7r/05) 


X cos (f sin(7r/a) — 


7r(^-l) 


a 


n—i 


(1 -n)/Q 


to {ccn-P-n + 2)Y 


where 




Mn-fi 


an-fl 


[s^ + iy 27ZijBra„ (i“ + l) 


-ds = ,^,{K.Y^ 


and where ci is given by ci = (—1)^L ji+...+j„=l ^1, ’ ■ ■ ■' 




b*j = 5o.j + q ■'(l-5oj) 


(!) 

q 

0 

0 ... 0 

(!) 

(!) 

q 

0 ... 0 

(!) 

(/i) 

(Al) 

. q 

(^) 


(4.) 

.(!) 

L+i) 

Y> 


(73) 


(74) 


(75) 


7GNU{0},^>0. 
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The function ^ j (f) is completely monotone whenever a € (0,1], 0 < ay < 
j3 < 1 HaED , and therefore by the Bernstein Theorem l26ll the spectral function ^ (r) 
is non-negative for the same range of the parameters. 

Moreover, from the above formulae ( |7^ and ( |74| ), we also derive the following 
result. 

Theorem 4.3. We have that 


poo 

L 


r) dr = 


a e (0,1], 7> 0, 


«e(i,2],7 = «eN. 

Proof. We let f —> 0+ in the formulae ( |72| and ( |74l i obtaining 

poo 

«e(o,i], 7 >o. 


(76) 


(77) 


and 


1 = 




From this, the claim easily follows. 

Corollary 3. If y=\, formulae Czl and (Izl), reduce to 


ae (1,2], 7 = «eN. 

(78) 

□ 


[ Kair)dr 

Jo 


1, as (0,1], 
1-i, as(l,2]. 


(79) 


The kernel Ka (r) has been thoroughly studied in the literature (for more in-depth 
information see e.g. Gorenflo and Mainardi lfT3]| and the references therein) while the 
general spectral function p{r) has been recently extensively analyzed in Mainardi 
and Garrappa ED. 

We conclude by emphasizing that, if a S (0,1], 0 < ay < 1, r > 0, the kernel 


<iW = 


r“y-i (yarctan ) + ?r(l - ay)) 

71 (r^®-|-2r“cos(7i:a)-b 1)^^^ 


(80) 


is the density of a probability measure concentrated on the positive real line (see Figures 

[Dl^andigi. 


Acknowledgements 

Zivorad Tomovski is supported by the European Commission and the Croatian 
Ministry of Science, Education and Sports Co-Einancing Agreement No. 291823. In 
particular, Zivorad Tomovski acknowledges the Marie Curie project FP7-PEOPLE-2011- 
COFUND program NEWFELPRO Grant Agreement No. 37 - Anomalous diffusion. 


19 








l(r) |(r) 



Figure 1: Plot (left) and log-log-plot (right) of the function ( [80| ) for a = 0.4 and 
7 = (0.25,0.5,0.75,1,1.25). 



Figure 2: Plot (left) and log-log-plot (right) of the function ( |80| ) for 7 = 4 and a = 
(0.05,0.1,0.15,0.2,0.25). 




Figure 3: Plot (left) and log-log-plot (right) of the function ( [80l l for 7 = 0.8 and 
a = (0.5,0.6,0.7,0.8,0.9). 
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